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Abstract 

We present a method for the calculation of the 03/2 heat kernel co- 
efficient of the heat operator trace for a partial differential operator of 
Laplace type on a compact Riemannian manifold with oblique bound- 
ary conditions. Using special case evaluations, restrictions are put on 
the general form of the coefficients, which, supplemented by conformal 
transformation techniques, allows the entire smeared coefficient to be 
determined. 
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1 Introduction 



The general topic of heat-kernel expansions and eigenvalue asymptotics of 
an operator L on a D-dimensional Riemannian manifold Ai has been an 
important issue for more than 20 years. In mathematics this interest stems, 
in particular, from the well-known connection that exists between the heat- 
equation and the Atiyah-Singer index theorem [Q]. In physics the expansion 
stands out in different domains of quantum field theory, as it contains, for 
example, information on the scaling and divergence behaviour, [0, 0, ^, ^, 

Of particular interest is the case where L is a Laplacian-like operator. 
When the manifold Ai has a boundary dA4 one has to impose boundary 
conditions that guarantee the self-adjointness and the ellipticity of the oper- 
ator L. The traditional, and simplest, conditions are Dirichlet and Neumann, 
and also a generalization of the latter referred to as Robin conditions. In this 
case, the heat-kernel coefficients have both a volume and a boundary part 
[0, |]. Writing the expansion in its usual form 



Kit)r^ Y: a,t'~^/' (1) 



one has the split 



fc=0,l/2,l, 



ak = dx bk + dy Ck- (2) 

JM JdM 



Here, the volume part, 6^, does not depend on the boundary conditions, 
whereas Ck exhibits a nontrivial dependence on the chosen boundary condi- 
tion. 

Whereas the calculation of the volume part nowadays is nearly automatic 
[y, [ly, |1T|, the analysis of Ck is in general much more difficult. Only quite 
recently has the coefficient C2 for Dirichlet and for Robin boundary conditions 
been found |T^, 0, 0, 0, 0, 0. C5/2 for manifolds with totally geodesic 



boundaries is given in |jT9[ 



One of the approaches is based on the functorial methods most systemati- 
cally used by Branson and Gilkey |jl2[ . Conformal transformation techniques 



give relations between the numerical multipliers in the heat-kernel coeffi- 
cients. However, on its own, this method is unable to determine the coeffi- 
cients fully. Additional information is needed, coming from other functorial 



relations or special case calculations |12 
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In contrast to this traditional boundary conditions relatively little is 
known about oblique boundary conditions. These more general conditions 
takes the form 

b = Vn + i(r"v„ + v.r") - s, (3) 

and involves tangential (covariant) derivatives, V^, computed from the in- 
duced metric on the boundary, is a bundle endomorphism valued boundary 
vector field and S" is a hermitian bundle automorphism. Finally, Vat is the 
outward normal derivative at the boundary. In order to ensure symmetry of 
the operator L together with the boundary condition 

BV\,^ = 0, (4) 

on a section of some vector bundle, one has to impose (F")^ = — F" and 
= S. This kind of boundary conditions arises naturally if one requires 
invariance of the boundary conditions under infinitesimal diffeomorphisms 
21| , or Becchi-Rouet-Stora-Tyutin transformations Furthermore 



@ is suggested by self-adjointness theory [0, and string theory p6| 

Although these boundary conditions have been subject of classical analy- 
sis (see e.g. ^ ^) the explicit determination of the heat-kernel coef- 
ficients has hardly begun. Some of the lower coefficients (up to Oi) have been 



evaluated by McAvity and Osborn [24] using their extension of the recursion 



method developed by De Witt for closed Riemannian manifolds. Another 



approach has been expounded in M^ based on the functorial methods used 



in [jT2| . In [^] the coefficients in the very special case of a flat ambient man- 
ifold with a totally geodesic, flat boundary have been computed. However, 
the combined knowledge of the coefficient with a totally geodesic, flat bound- 
ary and the functorial method relations does not allow for the determination 
of the full (under certain assumptions, see section II) 03/2 coefficient. The 
reason behind this is that no knowledge of the extrinsic curvature terms at 
the boundary is obtained via this example. This hinders the determination 
of all the universal constants. 

It is here that the approach of special case evaluation on manifolds with 
non- vanishing extrinsic curvature becomes important. A specific class of 
manifolds with this property is the generalized cone, a particular curved 
manifold whose boundary is not geodesically embedded. For this kind of 
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manifold techniques have been developed earlier for the evaluation of heat- 



kernel coefficients and functional determinants (see also [0, ^ ^ ^ 
P7| ) and have been applied to Dirichlet and Robin boundary conditions. 
Recently the generalisation to smeared heat-kernel coefficients was found 



This generalisation turns out to be important because functorial tech- 
niques (apart from other things) yield relations between the smeared and 
non-smeared case. The information one can get on the "smeared side" (these 
are terms containing normal derivatives of the smearing function /) is crucial 
to find the full "non-smeared" side (even present for / = 1). This has been 
demonstrated by determining the full 05/2 coefficient containing the whole 
group of extrinsic curvature terms PSI, ffT|. 



When employing the formalism of [32, 33 1 to oblique boundary conditions 



the tangential derivatives in (|^) cause added complication. For the case of 
the 4 dimensional ball these have been overcome in |]39| and the coefficients 



up to a2 were found for constant F* (for new considerations on this specific 
manifold see |^0[)- Here we will continue this analysis by generalizing it to 



the arbitrary dimensional as well as to the smeared case. As a result enough 
universal constants will emerge from special cases in order to find the full 
O3/2 coefficient for oblique boundary conditions. 

The organisation of the paper is as follows. In the next section we state 
in detail under which assumptions we are going to determine the 03/2 coef- 
ficient (purely Abelian problem, covariantly constant F*). For this situation 
the general form of the coefficients has been stated in and we explain 
which universal constants and which relations among them can be obtained 
from the generalized cone. In section 3 the explicit calculation on the cone 
is performed and the information predicted found. Having this information 
at hand the functorial techniques are applied to this boundary condi- 



tion To obtain as much information as possible even for the case of 

covariantly constant F* the structure and conformal properties of the non- 
covariantly constant F* will be needed and displayed. At the end, we will see 
that one relation is missing which will be obtained in section 5 by dealing 
with the manifold x T'^~^ {B"^ being the two dimensional ball and T*^"^ 
the {d — l)-dimensional torus). The conclusions summarize the main ideas 
and give further possible applications. 
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2 Restrictions from the generalized cone 



Before actually doing the calculation on the generalized cone we discuss the 
general structure of the coefficients and see what restrictions can be 
obtained from the coefficients found on our specific manifold. First one has 
to state clearly the assumptions under which the structure of the coefficients 
is formulated. We follow here |]31]], otherwise the situation is considerably 
more complicated [42|. The assumptions are as follows: 

(i) The problem is purely Abelian, i.e. the matrices F* commute : [F*, F-^] = 

0. 

(ii) The matrix F^ = /ij^FT-^ which automatically commutes with F-^ by 
virtue of (i), commutes also with the matrix S: [F^, S] = 0. 

(iii) The matrices F* are covariantly constant with respect to the (induced) 
connection on the boundary: VjF-' = 0. 

Under these assumptions we consider the Laplace-like operator 



L = -g'^ViVj - E 



(5) 



where E is an endomorphism of the smooth vector bundle V over A4 and V 
is a connection, together with the boundary conditions (^). Then the general 
form of the heat-kernel coefficients is 



ai/2(/) 



(47r 



|-(D-l)/2 

.1 



Tri5f)[dM] 



(47r)-^/2-{Tr(6/E + /i?)[A^] 
o 

+Tr {fiboK + b2S) + hif,N + faiKabV'T'] [dM\] 



(6) 
(7) 



03/2 (/) = (47r) 



-(D-l)/2_ 



-Tr 



ab 



384 

+C7SK + csS^) + f-wic^K + CgS) + Cef-NN 
+Tr 



[dM] 



(8) 



Here and in the following /[TW] = Jj^dxf{x) and f[dM.] = jQj^dyf{y), 
with dx and dy being the Riemannian volume elements of Ai and dAi. In 
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addition, the semi-colon denotes differentiation with respect to the Levi- 
Civita connection of A^, Q is the connection of V and Kab the extrinsic 
curvature. Finally, our sign convention for the Riemann tensor is R^jki = 
-^k,i + ^i,k + Kk^li - Ki^]k (see for example H). 



we 



Although ai/2 and ai have been determined previously [|^, 0, |39 
have included them in the list to explain clearly our procedure. The terms in 
a3/2 are grouped together such that the first two lines, Cq up to Cg, contain the 
type of geometric invariants already present for Robin boundary conditions, 
whereas all the other terms are due only to the tangential derivatives in the 
boundary condition. 

Our aim in the next section will be to put restrictions on the universal 
constants of eqs. (§P-(§P by calculating the coefficients of the conformal 
Laplacian {E = —{d — l)i?/(4(i)) on the bounded generalized cone. By 
this we mean the D = {d + 1) -dimensional space M. = I ^ M with the 
hyperspherical metric |Q 

ds^ = dr^ + r^dT.^, (9) 

where rfS^ is the metric on the manifold Af and r runs from to 1. A/" will 
be referred to as the base of the cone. If it has no boundary then it is the 
boundary of Ai. 

It is clear that a special case calculation will be simplified considerably by 
taking a constant F", say F*^ = ig, with the real constant g. In order that this 
is covariantly constant one might think of taking a fiat base manifold A/". The 
most natural such manifold where much is known about all required spectral 
propoerties is the torus. We thus choose A/" = T'^, namely the equilateral 
d-dimensional torus with perimeter L = 2tt and metric (iS^ = dxf -|- ... -|- dx'^. 
Its volume is vol(T'^) = (2ttY and the basic geometrical tensors read 

Furthermore, we choose a specific smearing function / = /(r) which will 
allow the calculation to be effected but which contains nevertheless all the 
information concerning the universal constants one can obtain. A possible 
choice is |Q (see section 4 for the applications in the given context) 

f{r) = h + hr' + f2r\ 
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For this special setting, the coefficients will have the following appearance: 
(47r)'^/2 



(47r)^/^ 
(47r)'^/2 



«l/2(/) 

6ai(/) 



5/(1) 



fo 



fi 



+ 



f2 



(27r)^ 



384a3/2(/) 



d-1 d+l d + 3 
+bofil)d + hMl) + 62/(1)5 - aJil)g 

/(I) \co{d - 1)74 - cid{d - 1) + carf^ 



(10) 
(11) 



-Ai^fif^ - /i2fi'7l - d) 



C4(i + cySci 



(12) 



+ /;Af(l) Csd + CgS* - /3i5f^ + Ce/; 



Thus, by doing the calculation on the manifold A4 = I x T'^ with /(r) = 
fo + /i^^ + /2^^ and by comparing terms containing a specific number of 
normal derivatives of / together with a fixed number of powers in d and 5" 
the following information can be extracted, 

ai/2 6, 
ai bo,bi,b2,ai, 

^3/2 C3 - Ci + Co/4, r^((J3 - 122) +C4 + C1- Co/2,C5,C6,C7,C8, Cg,/?!, Ai, 

rV2 + rV4 + rV2 + co/4. 

The considerable amount of information that one derives from this example 
is apparent, ai/2 and ai are completely determined without any additional 
input and from 03/2 one gets 10 of 18 unknowns. There is good hope that 
the remaining information can be found from functorial techniques. 

Having, therefore, a good motivation, we can embark on special case 
calculations and see afterwards if the functorial relations can complete the 
information. 



3 Non-smeared generalised cone calculation 

So let us turn to the spectral analysis of the conformal Laplacian on the 
described generalized cone together with the boundary condition (|^). The 
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conformal Laplacian is 

l- 

Ad dr"^ ' r dr ' Ar'^ ' r 



92 d d (d-lf 1 , 

- -^R = —,+z-. + + -2^^ (13) 



with eigenf unctions 



^(d-i)/2 exp{i(xini + .■. + x^nrf)}, n^TL. (14) 
The index z/ equals 

^=H + -+^9'^' (15) 
and the eigenvalues a are determined through (^ by 

aJl{a) + {u + gnd)Ju{(y) = 0. (16) 

Here u = l- D/2 - S. 

For the determination of the heat-kernel coefficients we follow the ap- 
proach developed in The basic object is the zeta function of Ai, 

c(^) = E«~'^ (17) 

and the relation 

ak/2 = r((D - k)/2)ResCM{{D - k)/2) (18) 

between the coefficients and the zeta function is used. In addition the Epstein 
type zeta function defined by 

n£ZL /{O} 

will turn out to be very useful. Obviously it is connected with the spec- 
trum of the Laplacian on the base manifold A/", the n^-powers arise from the 
tangential derivatives in (^. 

The starting point of the analysis of Cm is the contour integral represen- 
tation, 

Cis)= E [ ^k-'^^\nikJlik) + iu + gn,)Mk)), (20) 
Jy 2m dk 
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where 7 must enclose all the solutions of (|T6|) on the positive real axis. It is 
the appearance of the in the last term that causes the added complications 
compared to Q . 

In the following analysis the index z/ = would require a separate treat- 
ment. Its contribution has the rightmost pole at s = 1/2 because it is 
associated with the zeta function of a second order differential operator in 
one dimension. Because we are dealing with arbitrary dimensions this pole 
(and all other poles to the left of it) are irrelevant for our goal due to the 
relation (|18|) . For convenience therefore, we will continue without including 
the 1^ = contribution and will still use the same notation, C{s). 

Shifting the countour to the imaginary axis, the zeta function (with the 
zero mode i/ = omitted, as explained) reads 



C(^) 



sm ITS 



E 



n€2L'^/{d} 



dz 



zv] 



log2 

oz 



zvll{zv) + {u + gnd)Iu 



zv] 



As shown in detail in [^21, |33[, the heat-kernel coefficients are determined 
solely by the asymptotic contributions of the Bessel functions as — *• cx3. In 
the given consideration more care is needed since terms like Ud/v have to be 
counted as of order z/°. 



Using the uniform asymptotic expansion of the Bessel function 
encounters the expression 



one 



In 1 + 



1 + 



1 p 00 

E 

k=l 



Vk{t) 



1^ V V ) Z^^ \) 



00 

E 



Tj{u,g,t) 



f22) 



whereby the Tj are defined and t 
Uk and Vk, see . 

Asymptotically one finds 



1 / ^/l + z^. For the Olver polynomials, 



Cis) = A.^is) + Ao{s) + A4s) + J2A^ 



(23) 



(21) 



where A_i{s) 
namely 



and A 



A-iis] 



are the same as in Robin boundary conditions 
1 r(s-l/2) 



4^^ r(s + l) 



-Eo{s~ 1/2), 



(24) 
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Ms) = Ieo{s). (25) 
The new quantities are 

and 

^.W— E f^. (27) 

In order to proceed it is convenient to express Tj as the finite sum 

= /a Ac 777^776' (28) 

with 5 = grid/v. The /f"'^^^ are easily determined via an algebraic computer 
programme. 

The next steps are to perform the 2;-integrations by the identity, 

^ (i + 5t)^"2 r(y) ,^0^ k\r{{x + k)/2) ''^''^ 

and then do the n-summation to write everything in terms of the Epstein 
functions ([T9|). Performing these steps one gets first 

For c even and c odd in (|28| ) slightly different representations appear such 
that some more notation is unfortunately necessary. We write 

a,b,c 

and get 



^.(^) = E/iK'^(^) (31) 



_ 1 ^ T{b + 2n) Tis + a/2 + n 
^3 - ~TT7{ 2^ 



T{s)^,V{b)V{2n + l) V{a/2 + n) 

^2n+c(s + n + (j + c)/2)/'^+^ (32) 
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for c even and 



^ 1 ^ r(& + 2n + l)r(5 + (a + l)/2 + n) 
' r(.)„^,r(6)(2n + l)! r((a+l)/2 + n) 

E2n+c+i{s + n + {j + c+ l)/2)/"+^+\ (33) 

for c odd. 

We need the residues of v4_i, Aq, and A^, but this is not too difficult, 
because the Epstein zeta functions are very well studied objects. For us the 
relevant properties are 

En{s) = for n odd (34) 
and for s = I + d/2 the residue is 

(d-l)/2p/; I I /2^ 

A nice feature of the calculation is, that when using the above results (Q) 
and (|35D in {^) and (|3l|). Res A+{{D - k)/2) and Res Aj{{D - k)/2) reduce 
to the series representation of the generalized hypergeometric function 



{ai)k{a2)k---{ap)k 
to {PMP2)k..mk kV 



pFq(ai, 02, •-, ap] A, /?2, Pq] z) ^ 77TY77TY— t^^Vtt- (^6) 



For example, for only one contributions arises which, usefully normalized, 
reads 

r((D-l)/2)L_L-Res A4(D-l)/2) 
= ^ {2^1(1/2, rf/2;rf/2;/)-l} 

= ^{(1-/^/^-1}- (37) 

In this case the intermediate step in terms of the hypergeometric function is 
artificial of course but useful in general. 

In order to give the contributions of A'^'^''^, we have to distinguish between 
even and odd j. Although a bit lengthy we find it useful to state these results. 
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Let us stress, that we have already, in principle, determined an arbitrary 
number of heat-kernel coefficients for the Laplacian on a generalized cone 
with oblique boundary boundary conditions in the non-smeared case. For c 
odd and j odd we find 

r{{D - 1 - j)/2)^|^Res Af^%D - 1 - j)/2) = 

g r(i + ./2) 

r((a + l)/2)r(6) 

6-1 



^ + c d . , 



a—J—c 
2 



.M^J 9^ (38) 

3F2(1, {d + a + 1- j)/2, 1 + c/2; (a + l)/2, (L> + c)/2; /), 

which contributes to the coefficient a(i+j)/2. For the specific values of 6, c,j 
and A; needed the hypergeometric function always reduces to a simple alge- 
braic or a hyperbolic function. The above result neatly summarizes all this 
information in one equation. 

For c odd and j even A"'*''^ also contributes to the coefficient a(j+i)/2 and 
the relevant result is 1/ {2^/71) times the above (note that in this case the 
normalization is (477)*^/^). 

Furthermore, for c even and j odd the analogous result is 



(47r) 



D/2 



T{{D - 1 - j)/2)^^^Res Af'^{{D - l-j)/2) 
r((l + c)/2) 



r(a/2)r(6) 



6-1 



2 

3F2(1, (d + a- j)/2, (1 + c)/2; a/2, (d + c)/2; /), 

For c even and j even the same rules as above hold, furthermore the same 
comments. 

The above results allow for a direct evaluation of the coefficients by an 
algebraic computer program such as Mathematica. However, before stating 
the results let us describe the necessary modifications when dealing with the 
smeared case. 
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4 Generalization to the smeared heat-kernel 
coefficients 

The inclusion of a smearing function that depends only on the radial variable 
results in the smeared zeta function, 

dk 



2m 



-k 



-2s 



' dr f{r)Jl{kr)r ^ ln{kJl{k) + {u + gn,) J,{k)) .{40) 
ok 



af;s)= E 

ne2L'^/{d} 

For /(r) a polynomial, we have shown how to analyse ( ^OD for Dirichlet and 
Robin boundary conditions in and so for the general procedure see this 
reference. The generalisation to oblique boundary conditions is obtained 
here. 

The bar in (^OD signifies normalized. Explicitly 

V2k 



Mkr) 



{{u + gndY + k^-p^Y/^J,{k) 



Mkr) 



For 



N 



2n 



n=0 



we need normalization integrals of the type 

S\l + 2p] = dr Jl{ar)r^^^P . 
Jo 



(41) 



(42) 



(43) 



These can be treated using Schafheitlin's reduction formula which for 
the present case gives the recursion 

2p — 



S[l + 2p\ 



S[2p - 1] 

2p{u + p) 

2p + 1 ' + {u + gud)'^ 



2p + l 
1 



1 + 



2 ' 



(44) 



starting with ^[l] = 1. So ^[l + 2p\ has the following form 

V / ,,\2m ra 



m=0 



V 



1=0 
1 



a"^ + {u + grid)'^ — 



m=0 



a 



2m n 



-21 



ml 



(45) 



1=0 
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The numerical coefficients 7^; and are easily determined recursively. As 
a result, apart from characteristic differences, the smeared zeta function 
s) takes a similar form as Cm{^)- 
It will turn out convenient to devide CmUiS) into different pieces char- 
acterized below. First, respecting the structure in (|45|), we define 

dk 



= EE7L E z." / 

-ln{kJl{k) + {u + gnd)Ju{k)) (46) 



d_ 



m,s) = EEC. E 

n 7 r. 7 -'■7 



27rz (fc2 + (m + gndY — 



d 

—ln{kJl{k) + {u + gnd)Ju{k)), (47) 

where the contour 7 has to be chosen so as to enclose the zeros of only 

kJl{k) + {u + gnd)Jy{k). Thus the poles of ^[l + 2p], located at = 



±y — (u + gn^y , must be outside the contour. It is important to lo- 
cate the contour properly because, when deforming it to the imaginary axis. 



contributions from the pole at = yz^^ — (u + gudY arise. 

The index p referes to the fact that these are the contributions coming 
from the power r^^ in (^). In order to obtain the full zeta function, the 
J2p=o fpC^ has to be done. 

The ffist piece, may be given just by inspection. Comparing (^) 
with the non-smeared zeta function ( pO| ) the contour integral is the same 
as previously once s ^ s + m has been put. Due to the additional factor 
^2m-2i argument of the base zeta function has to be raised by / — m. For 
explanatory purposes let us give as an explicit example 

In exactly the same way, y4o(/, s), A']_{f,s) and y4j(/, s) are obtained from 
(|25D, (|30D, (|32D and (^). This is the stage where the properties ( p4D and 
(|35|) are used and the contributions to the heat-kernel coefficients in terms 
of hypergeometric function emerge. They will not be displayed, however. 
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explicitly, because the structure is the one already seen in ( pSD and the way 
they are obtained is identical to the procedure described in section 3. 

We continue with the analysis of where several additional complica- 
tions occur. Shifting the contour to the imaginary axis we get the pieces 

p—l m 

m,s) = !^ ^ EEC(-ir^-^^-^^ (49) 

i-oo ^-2s-2m Q 

/ dz - — ■ 2N ^^^ ((^ + 9nd)Iu{zv) + zvI'A^z)) 

Jo [u + gud) — I' [I + z^) oz 

-| p— 1 m 

CLM.s) = -^EEC. E i^'-'^i^'-iu + gudrr (50) 
d 

— In {kJl{k) + {u + gnd)Uk))\k=^ ,.^^^+g^^)2 . 



the last one arising on moving the contour over the pole aX. k = ^ v"^ — [u + gudY- 

In dealing with (^{f,s) one can use, as done after eq. (^), the uni- 
form asymptotics of the Bessel functions and define analogously to (|23| ) the 
asymptotic contributions Ai^sif, s). We will illustrate the calculation by deal- 
ing with 



A\sif,s) = !^ ^ EEc.(-i; 



p—l m 



oo ^—2s—2m—l 



dz - (1 + ^2)1/2 

{u + gudY -u^{l + z^y ' ■ 



By using the expansion 



1 _g {u + gudf' 



[U + gUdf - 1/2(1 + Z^) Z/2i+2(l + ^2)i+l 



the above integrals are recognised as representations of Beta functions 
As an intermediate result one gets 



A\,{f,s) = EE EC 



°° T{s + i + m + l/2) 

.=0^=0 /=o r(. + m + l)r(^ + l/2) 



2-^ ,,2s+2i+2i+l ■ I 
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Complicated as the expression (^TJ) is, we have to remind ourselves of our 
initial goal, namely, the determination of the heat-kernel coefficients up to 
03/2- Thus we need to determine only the residues of (^) at s = -D/2, {D — 
l)/2,D/2- 1 and {D-3)/2. 

In general, the numbers contain terms independent of and linear 
in Ud, 

Furthermore it is clear that the higher the power of the more to the 
right the pole of the associated term will be. Thus in addition consider the 
expansion 

{u + gn^r = g^'nf + 2m/- + O(nr^) 

in powers of ria. With eq. (^) it is then obvious that the rightmost pole 
in A^^^g{f,s) due to the term is situated at s = (D — 4)/2 and 

contributes only to A2. For our immediate purposes it is thus sufficient to 
take into consideration only the above two terms. As a result 



with 



A\sif,s) = Ffif,s) + F!if,s) 



_ l^^vV^^ r(. + ^ + m + l/2) 
' ^ 2 to h h '"'^ r(. + m + l)r(z + 1/2) 

g''E2i{s + 1 + 1 + 1/2), (52) 



'^U^oto r(. + m + 1)1(^ + 1/2) 

g'''E2i{s + 1 + 1 + 1/2). (53) 

Use of the residues of the base zeta function then easily gives the following 
normalized contributions (we use the notation = for / < 0) 



r{D/2-k)^—L^ResF[{f,D/2-k)= E 



m+l— fe 



(2vr)^ ^ 1 V. > - ^ " ^fr^i {D/2 - k)^^, 

2Fi{l,D/2- k + m + 1/2, d/2; g^) (54) 
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T{D/2 - k)^-^^Res D/2 - k) 

dg 



u 



E 



=fc-i 



l,L)/2-fc + m + l/2,rf/2;^^ 



cp {d/'^)m+l-k 

-{fc-l)lp/2_A;)„,+i 
(55) 



For our purposes, only A; = 1 is relevant, but we have given these general 
results to show that in principle one could go further. 

The same kind of argument allows one to show that the relevant parts 
in the other Af^f^s) can all be representated in terms of hypergeometric 
functions. 

Finally we are left with the treatment of Cshiftif^^)^ ^l- (HO)- Here, 
instead of the uniform asymptotics of I,y we need it for the function J^. The 
45| (see also ||38|) suggests for z/ — oo. 



expansion given m 
d 



dk 



ln{kJl{k) + {u + gnd)Ju{k)) 



1=0 



2«+l 



(56) 



with the coefficients ti to be determined. Here the problem appears, that 
every value of / contributes to the pole of Cshiftify ^) already at s = (D — 1)/2. 
Thus, the asymptotic expansion of the left hand side to any power in ria/i' is 
needed, apparently an extremely difficult problem on asymptotics of special 
functions. However, given that Cshiftif^ ^) contributes only for p > we will 
show in the next section how to circumvent a direct evaluation of eq. (|56D. 



5 Results and conformal techniques 

After having shown in detail several aspects of the special case calculation 
on the generalized cone let us collect the information about the universal 



constants appearing in eqs. ([lOD - (0). 

In order to answer the open question about eq. (|5^) let us consider first 
ai/2 and take /(r) = /q. For this case only Ao{s) and A+(s) contributes and 
the answer is 

. f,l (_^ _ i) . (57) 

Comparison with (|r^) gives the correct universal constant 

^=U > ^ - (58) 
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and the present special case evaluation determines the ai/2 coefficient for a 
general manifold (which was clear of course). Thus, for the case of the cone, 
taking /(r) = /o + fir"^ we know that 

\^av.(/) = (/o + /Oj(^-l]. 

Taking into account all terms but the unknown contribution from CsMftif^ ^) 
we find 

1/2 \ 

Since the ffist piece is the correct answer, the last piece has to be cancelled 
by the contribution of Cshiftify Use of the expansion (^) and (|50D gives 

V 2 / i-^-l 

which leads to the condition 

°° °° (1) (¥) 1 



j=0 «=0 , V ■ / • , , 



In order that the left hand side should be nothing other than a complicated 
series expansion of the right hand side one has to conclude that ti = (1/2)^//! 
and one finds the surprisingly easy expansion 



d 

— ln[U'^[k) + [u + gn,)) ^^^/-.-^-^ _ ^ — 

But having expansion (|59D at hand the contribution of CsM ftify ^) any pole 
can be determined so that now for the cone complete knowledge for ai and 
03/2 is available. 

We turn now to ai. Dealing ffist with /(r) = /o, the linear term in d 
defines bo, the linear term in 5* defines 62? the term independent of d and 
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S defines ai. Dealing afterwards with /(r) = /ir^ the additional piece is 
immediately identified with 61. As a result we obtain the correct answer 



60 = 2-6 

61 = 3- 



1 



1 + r^ 



+ 



7^ , 



6 ArcTanh(\/— r^) 



12 



l + P' 

_6_ / 1 

f2 l~i + r2 



+ 



ArcTanh(y^ 

7^ 



which shows that the ideas involved in our special case calculation are in- 
deed correct. Especially the example once more determines the complete 
coefficient for a general smooth manifold with boundary. 

Proceeding in the same way as described for ai, we obtain the following 
universal constants for 03/2, 



C5 



1 



C6 = 



8 (32 



144 



^0 



160 



32 



-1 + 



1 



1 



32 



(60) 



-36 + 



52 



32 (5 - 8 (-1 + v^TTp) - 4 (-4 + 3 VTTT^] 



cr = 



C8 = 



192 (1 - v^TTp - (-2 + y/T+T^ 



r2(i + r2)5 



192 



:i + r2)i' 



-192 



Pi 



-32 (5r4-8 (-1 + Vi + r2) -4r2 (-4 + 3^rTr2 



(61) 

(62) 
(63) 
(64) 

(65) 
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192 



Ai 

r4(i + r2)2 

respectively the following relations among them 

C3-C1+ Co/4 = p4(x + p)3/2 

-32 (-1 + Vl + T^) - 5 (-16 + 3 Vi + r 



(66) 



(240 - 224 Vl + P) + (336 - 207 Vl + P 

(67) 



r^((T3 - /ia) + C4 + Ci - Co/2 



32 (-1 + Vl + P 



r4(i + r2)3/2 
(-48 + 7 Vi + p) + 16 (-10 + 9 VTTT^ 
(-192 + 119 Vi + P) 



(68) 



1 \ _ 8 



32 (-1 + Vi + r^ 



f-32 + 3 Vi + tA + 8 f-15 + 13 Vi + r2 



+3r^ (-42 + 25Vr + T2 



(69) 



Where applicable, in the limit F ^ the results for Robin boundary condi- 
tions are reproduced check. 



This serves as a very good input for applying the techniques of [|T2[ . First 
we use a result on product manifolds [O, which in our case gives 



Co 



Cl 



96 
16 



-1 



-1 + 



2 



Together with eq. (0) this also determines C3, 



(70) 
(71) 



C3 



F^ (240 - 224 Vl + F2) + F^ (320 - 199 VlTf' 



F4(l + F2)3/2 
+F^ (64 - 7 v^ITn) - 32 (-1 + Vl + F^ 



(72) 
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The remaining task is to apply the functorial techniques of |T^. For de- 
tails of the technique itself see this reference and for the modifications when 
tangential derivatives are involved see pi . 



The basic equations are the conformal-variation formulae 

d 
de 



xl-o«n/2(l,e"2^^L) = {D-n)an/2{f,L), (73) 



^^U=oan/2{e-^'^H,e-^'^L) = for D = n + 2, (74) 

with an arbitrary smooth function H. For a collection of variational formulae 
again see [Q. The additional relation r*(e) = e'^-^F* is given in ||3l|. Setting 
to zero the coefficients of all terms in ( [73| ) for n = 3 one obtains, for example 
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Term Coefficient 

f;NN = i(D - 2)C0 - 2{D - 1)C1 -{D- 1)C2 -{D- 3)C6 - FVi 

Kf,N = i(D - 2)co - 2{D - l)ci - C2 + 2{D - l)c3 + 2c4 
- liD - 2)cj -{D- 3)C5 + FV3 - FV2 

The first of these determines C2 and /ii, namely 

96 f 2 + F2 - 2 VT+T^) 
/xi = ^ , ^. (76) 

The second together with (|^) gives 

= 2 / , / _ U - ^) + 32 f-1 + (77) 



In addition we get 



h = 0. (78) 



Disappointing as it is, under the given assumptions these are the only new 
universal constants the functorial techniques yield. But due to the restrictions 
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imposed we have not yet exploited all information available. For example one 
has the variational formula 

' ,R\ipV = -2fR\i^rT^ -{D-?,)K,f^Vf,r, (79) 



de 



Up to now we have not compared coefficients involving tangential deriva- 
tives, because for covariantly constant F* these pieces integrate to zero. But 
if we relax the condition of covariantly constant F*, eq. (|7^) shows that two 
additional eqs. for the universal constant ^2 would arise. (In contrast the 
variations of the terms associated with the missing (T2, (T3 and a 4^ do not con- 
tain any tangential derivatives.) In order to exploit this observation we have 
to generalize eq. (||) to when VjF* 7^ 0. This is done as usual by building up 
all possible independent geometrical terms with certain homogeneity prop- 
erties fl^. For the boundary conditions under consideration eq. (|^) has to 
be supplemented by the following terms, 

^^a72(/) = Tr[/(7iF\/,^+72r^/^, + 73r^F^', 

+ 75r,/F^ + 76F^/,.fcF^F'= 

+7rF^/fc;,F^F'= + 78F^/,^^F,F'= + 79F^,F'^,F^F*^ 

7ior,,,F,F^F'^ + 7nr^/'^,F^F'=F,F0][aA^]. (80) 

The term F*. • F-' is not added because due to the Gauss-Codacci relation one 
■J * 

has 

so that this term is linearly dependent on the others already displayed. Be- 
cause of the simple conformal transformation property of F* it is relatively 
easy to find the variational formulas of all invariants in eq. ([80|). Even 
though we know none of the 7^, setting to zero the coefficients of the tangen- 
tial derivatives terms in ([79|), we find the unambigous answer 

^l2 = 0. (81) 

As a consequence eq. ( |60D shows 
1 



0-3 



re (1 + r2)3/2 



32 (-5 F'' + 8 (-1 + VT+T^) + 6 F^ (-5 + 3 VTTT^) 
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(-30 + 26 Vl + r2)) 



This is really all we can get from the lemmas and the specific example 
of the generalized cone, I x T'^. We can obtain information only about the 
combination r^a2 + (T4. Thus some additional input is needed to accomplish 
the goal of finding all universal constants in (H). 

6 Oblique boundary conditions on x T^''^ 

One possibility of finding the remaining information is to look for an example 
which is able to separate the contributions of (f^yTT-^)^ and KijKlT^VK 
The reason that all types of generalized cones with metric (|^) fail to do so, 
is that Kl = Sj. As a result, the contraction KijT^T^ equals and also 
KijKiV'V^ = T\ Having a term like a2iT^)T^ + a4{T^)T^ = g{T^), with 
g{T^) a known function of F^, there is no possibility uniquely determining 

or (T4 because, as indicated, these also depend on F^. It is clear that this 
problem is not a result of having chosen only one non-vanishing component 
F*. For a generalized cone, these properties are generic. So we are forced to 
leave this class of examples. 

The difference between the invariants associated with a2 and is that 
the first contains fourth powers of F* whereas the second one only squares. 
If we deal with two instead of one nonvanishing component, say F*^ = and 

= dd, i d, gd, g constants, and if the extrinsic curvature could be a 
projector on one of them, say K^d = 1, Kij = for (z,j) 7^ {d,d), then 
[KijT^r^)^ = g^ and KijKlV^T^ = g"^. By simply comparing powers of gd the 
universal constant 02 could be determined being the only one with g'^. 

Using the metric 

ds^ = dr^ + dT? 

and Kab = — Tab, Tab being the Christoffel symbol, it is seen that keeping the 
manifold topologically as J x T*^, the metric 

dT? = dx\ + ... + dx'^^_i + r^c/x^ 

will have the above property. This is clearly the flat manifold x T'^^^ and 
the eigenvalue problem is easily solved. With the notation n| = n\ + ...+v?^_i 
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the eigenf unctions are 



with eigenvalues o? . The boundary condition takes the form 

^d^ - nlJ\^J^\ja^ - nl) + {g^nd + gn - S)J\n^\{\ja^-nl) = 0, (83) 

where we have used n = rii (the result is the same for any i & {1, ...,d — 1}). 

Our main interest is to determine a2 and a 4^. The calculation involving 
two tangential derivatives will be seen to be sufficiently difficult so that we 
will restrict ourselves to what is strictly necessary, namely we will not bother 
to do the smeared calculation but will be content with the sufficient choice 
/(r) = 1. The information derived about C3 + 04,07, C8,Ai will serve as a 
further check of the previous calculation, the new quantities, (T2 and (T3 + (T4 
will complete our analysis of 03/2 (/), for covariantly constant F*. 

One can procede very much as before. Starting with a contour represen- 
tation similar to eq. (^) and shifting the contour to the imaginary axis, one 
gets 

as) = E / ,rfM^'-r^?)-^T^ln K„,|(fc) + [^7.n, + ^n-S]J|„,|(fc) .(84) 

It is seen, that flf acts effectively as a mass of the field and comparing with 
the previously treated example plays the role of u. 

The general procedure of dealing with ({s) is the same as in section 3. 
However several complications arise, and the situation is sufficiently different 
so as to warrant further describtion. 

One starts from the uniform asymptotic expansion of the Bessel function 
[H] and eq. (p2D is found with the replacements gn^ gd^d + gn, u — *• —S 



and z/ — > \nd\. Asymptotically one again finds eq. (^) with the characteristic 
differences already described. The several new features arising are now dealt 
with by looking at 

Ms) = ±11 E' r dz[z^nl-nlr (85) 

gn + gdUd 1 \ 
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The integral is nothing but a hypergeometric function |^ 



and we get 



A is) - 1 V( ,y r(. + (/ + l)/2) 

^ 2r(.)^^ r((/ + 3)/2) 

oo 

^ // + 3 /+1 / + 3 Ud 



(86) 



The apparent difficulty is to extract the meromorphic structure of multiple 
sums of hypergeometric functions. This is very effectively done by using the 
Mellin-Barnes integral representation of 2-^1 P6|| , 



r(7) 1 Tia + t)T{P + t)Ti-t)^ 



^^^^"'^'^'^^ = fRfM2^Loo'' fWTT) 

where the contour is chosen such that the poles of the function T{a + t) and 
T{P + t) lie to the left of the path of integration and the poles of the function 
r(— t) lie to the right of it. When using this integral representation it is seen 
that the sum over leads to (rf — l)-dimensional Epstein type zeta functions, 
whereas the sum over Ud gives a Riemann zeta function. The relevant zeta 
function is again of the Epstein type (|19D which we now write as 

EtMs)= E i^'tr'n'^- (87) 
As a result of the described steps one arrives at 

rioo 

— dtT{s + l + t)r(-t)C/j(-2t -21 + 2k)Et,2k{s + l + t), 

ZTTt J —ioo 

where the contour (depending on / and k) is such that the poles of (r lie to 
the right of the contour, the poles of Et^2k to the left of it. This Mellin-Barnes 
representation allows the meromorphic structure of A^{s) to be read off by 
closing the contour to the left. We then encounter poles of T{s + I + t) at 
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t = —s — I — m, m E INo with residues r(s + / + m)^ij(2s + 2k + 2m)Et^2k{jn)- 
The right-most pole hes at s = 1/2 and it is clear that the poles of the F- 
function are irrelevant for our purposes. However, the pole of the Epstein 
function is situated at t = {d—l)/2—l — s + k and gives relevant contributions. 
Keeping only these terms, 

vr('^-^^/'^l ' /2/ \r(fc+l/2)^^ , , ^^^^^ 

C«(2s - d + l)g''g'f- (89) 

The right most pole at s = d/2 comes from the Riemann zeta function. Using 
the relations |^ 

T(x) 



r(2x) 22^-ir(x + 1/2) 

(2Z)! _ (2/-1)! _ r(/ + l/2) 

~ 2^ ~ 77^ 



and 



one gets 

Res A^{d/2) = {(1 - / - g'.r'/' - 1} . (90) 

The calculation shows the manner in which = g'^ + g^ is built up. Together 
with the contribution of Aq{s) one finds the correct coefficient ai/2- 

The next pole in ( [74D at s = (c? — l)/2 comes from the F-function for 
k = l. Then 

Res A^{{d - l)/2) = -^r^l^^?^ {(1 - - 1} • (91) 

This piece is cancelled by another contribution, but the example shows the 
way other contributions than F^ appear. 

There are no further (interesting) poles due to the zeroes of Cr{s) at 
s = —2m, m G IN. 

The basic characteristics, namely that the integrals are hypergeometric 
functions and that the eigenvalue sums may be dealt with by Mellin-Barnes 
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integral representations of these, are present for all other Aj{s). The final 
results can all be written in terms of 2-^1 or 3F2 or their derivatives, very 
much as in section 3. For the sake of space it is impossible to give further 
details. 

In summary we confirm our results for the universal constants C3+C4, C7, Cg 
and Ai- Most important, we determine the remaining constants to be 



0-2 



r8(i + r2) 



-48 (-51^ 



16 



(-30 + 16 Vi + r 



4 + Vi + r2) + 8 (-5 + 4 Vi + r2 

(92) 



0-4 



32 (-r^ + 16 (-1 + VI + r2) + 2 r2 (^-7 + 3 VTTT^ 



re vi + P 

and have thereby achieved the goal of determining 03/2(7) in eq. (|^). 



(93) 



7 Conclusions 

In this article we have developed a technique for the calculation of smeared 
heat-kernel coefficients on generalized cones or manifolds of the type x 
rpD-n operators of Laplace type with oblique boundary conditions. These 
boundary conditions arise in response to questions in quantum gravity, gauge 



theory and string theory [^, |J, ^ . The asymptotic properties of 

the generalized boundary conditions encoded in the asymptotic heat-kernel 
expansions are considerably more involved than the corresponding ones for 
the traditional conditions. This article makes an attempt to provide a prac- 
tical approach for the calculation of coefficients to any order needed although 
we suspect that the involved analysis, complicated results and restrictions in- 
dicate that further analysis along the lines of this paper should not be lightly 
undertaken without strong motivation. 

The direct calculation of higher coefficients for general curved manifolds 
with arbitrary smooth boundaries becomes very difficult and impractical. In 
the approach promoted here (see also [^) this analysis is avoided. Based on 
the observation that functorial methods give relations between the numerical 



multipliers in the heat-kernel coefficients |T2| the remaining task is to find as 
many multipliers as needed by other means. A rich source of information is 
special case calculation. Done in a systematic fashion in arbitrary dimensions 
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and with a smearing function as general as necessary, we have seen that 
simple comparison with the general form of the coefficient yields the encoded 
information on the universal constants. Several checks on the constants found 
are provided by the conformal relations as well as by the different examples 
treated. A systematic feature of the approach is the use of algebraic computer 
programs. 

It is clear that the approach can also be applied to the calculation of a2, 
however, more than 100 terms are involved and a very real additional effort is 
necessary as well as a substantial motivation. Generalization to covariantly 
nonconstant F* is desirable and can be attacked by taking different base 
manifolds (eg. a sphere) or by considering simple dependences of on one 
of the tangential variables. Finally, it is hoped that also for the non-Abelian 
setting new information can be obtained in the spirit of the present work. 
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